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Abstract 

At the end of their birth process, neutron stars can be subject to a magnetorotational instability in which a conversion of kinetic 
energy of differential rotation into radiation and kinetic energies is expected to occur at the Alfven timescale of few ms. This birth 
energy conversion predicts the observed large velocity of neutron stars if during the evolving of this instability the periods are of 
few ms and the magnetic fields reach values of 10 16 G. 
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A major unsolved problem in astrophysics is to explain why 
neutron stars exhibit space velocities well above those of their 
progenitor stars (Anderson & Lyne , 1983; Lai et ah, 2001; Lai 
, 2004; Fryer , 2004; Chatterjee et ah, 2005; Hobbs et ah, 2005; 
Postnov & Yungelson , 2006). Most neutron stars have been 
measured to possess high space velocities in the range of 100- 
1000 km/s, while their progenitor stars have velocities of the or¬ 
der of 10-20 km/s (Lyne & Lorimer , 1994; Arzoumanian et ah, 
2002; Fryer et ah, 1998; Hansen & Phinney , 1997; Lorimer 
et ah, 1997; Cordes & Chernoff , 1998). It is generally accepted 
that neutron stars receive a substantial kick at birth, which pro¬ 
duce their observed space velocities. However, the physical ori¬ 
gin of this kick is unclear. Some of the proposed kick mecha¬ 
nisms require large initial magnetic fields in the magnetar range 
(10 14 -10 16 G) and asymmetric emission of neutrinos (Kusenko 
& Segre , 1996, 1997; Lai & Qian , 1998; Maruyana et ah, 
2011). Other mechanisms require a rapid initial rotation to 
produce substantial kicks (Spruit & Phinney , 1998; Khokhlov 
et ah, 1999; Sawai et ah, 2008). Hydrodynamical models have 
also been suggested, which are based on recoil due to asymmet¬ 
ric supernovae (Burrows & Hayes , 1996; Scheck et ah, 2004; 
Burrows et ah, 2007; Nordhaus et ah, 2012). Recent models 
rely on the existence of topological vector currents (Charbon- 
neau & Zhitnitsky , 2010). The current distribution of the ob¬ 
served neutron star velocities seems to be Maxwellian, which 
points to a common acceleration mechanism (Hansen & Phin¬ 
ney , 1997; Hobbs et ah, 2005). The elusive kick mechanism 
appears to be connected with the possibility that during some 
stage of their birth process, neutron stars reach magnetic fields 
typical of magnetars and periods typical of millisecond pulsars 
(Usov , 1992; Duncan & Thompson , 1992; Thompson , 1994). 

In this paper we suggest that the kick velocity of neutron stars 
may arise from a magnetorotational instability (MRI) produced 
at the end of their birth process. More specifically, a newly-born 
neutron star can be subject to a MRI evolving at the Alfven time 
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of a few ms in which an emission of radiation energy accompa¬ 
nied by a gain of kinetic energy of translation can be produced 
at the expense of a loss of kinetic energy of differential rotation. 
If during the evolving of this birth MRI the period is of a few 
ms and the magnetic field reaches values of 10 16 G then we show 
that the gain of kinetic energy of translation can predict the ob¬ 
served large kick velocity of neutron stars. Our suggestion is 
supported by studies showing that there is an amplification of 
the magnetic held and a transference of angular momentum dur¬ 
ing the evolving of MRI in newly-born neutron stars (Akiyama 
et al., 2003; Thompson et ah, 2005; Masada et ah, 2012). Simu¬ 
lations have shown that birth MRI can generate magnetic fields 
of the order of 10 16 - 10 17 G in several ms (Thompson et ah, 
2005; Siegel et ah, 2013). 

Newly-born neutron stars are assumed to be highly convec¬ 
tive and differentially rotating hot fluids (Stergioulas , 2003; 
Yamada & Sawai , 2004), which can be subject to MRI evolv¬ 
ing at Alfven timescales (Duez et ah, 2006). During this birth 
MRI the star can convert its kinetic energy of differential rota¬ 
tion into magnetic energy (Akiyama et ah, 2003; Spruit , 2008). 
It is then plausible to assume that the star can lose kinetic en¬ 
ergy of differential rotation via radiation (a fraction of the total 
magnetic energy). Conceivable, in this stage the star can also 
gain kinetic energy of translation and this is the basic assump¬ 
tion of the model proposed here. We assume that at the end of 
a birth MRI, neutron stars experience the energy conversion 


d(Mv~ 

rad + dt\~Y~ 


d ladlQr 
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where P ra( j is the instantaneous radiated power; Mv 2 /2 is the 
kinetic energy; ad/D 2 /2 is the kinetic energy of differential ro¬ 
tation; v is the space velocity; a<± is a dimensionless constant pa¬ 
rameter accounting for the differential rotation (Spruit , 2008); 
I = 2MR 2 /5 is the moment of inertia with M and R being the 
mass and radius; and O = In)P is the angular velocity with P 
being the period. In Eq. (1) we have assumed that there are not 
significant changes of the parameter ad during the short period 
occurring the birth MRI. But in general ad varies with time. 
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According to Eq. (1) an emission of radiation energy and an 
increase of kinetic energy of translation occur at the expense 
of a loss of kinetic energy of differential rotation. We note 
that an equation similar to Eq. (1) [without era and with P ral j] 
associated with the asymmetric radiation from an off-centered 
magnetic dipole is the basis of the “rocket model” proposed by 
Harrison & Tademaru (1975). Another equation also simi¬ 
lar to Eq. (1) [without a<± and with P ra(1 estimated by an ex¬ 
ponential field decay law] has been considered to study both 
the birth accelerations of neutron stars ( Heras , 2013) and the 
birth-ultra-fast-magnetic-field decay of neutron stars ( Heras , 
2012). It is pertinent to say that the idea that newly-bom neu¬ 
tron stars would lose their rotational energy catastrophically on 
a timescale of seconds or less was suggested by Usov (1992). 

The instantaneous radiated power P ra d in Eq. (1) must express 
the idea that a MRI is responsible for a rapid exponential growth 
of the magnetic field (Akiyama et al., 2003; Yamada & Sawai 
, 2004; Siegel et ah, 2013). In this sense. Spruit (2008) has 
pointed out that some form of MRI occurring during a differen¬ 
tial rotation in the final stages of the core collapse phase may 
produce an exponential growth of the magnetic field and that 
once formed, the magnetic field is in risk of decaying again by 
magnetic instabilities. Without considering changes of kinetic 
energy, it is well-known that abrupt changes of kinetic energy 
of rotation produce abrupt changes of energy of radiation. The 
characteristic time of rotation changes is similar to the charac¬ 
teristic time of radiation changes. Therefore if the birth MRI 
occurs at Alfven times of ms then the radiative energy must be 
emitted in these times and therefore the increase and decrease 
of magnetic fields producing such a radiative energy must oc¬ 
cur at these times. Accordingly, we assume here the existence 
of a birth MRI responsible for a rapid exponential growth of the 
magnetic field, followed by an equally rapid exponential decay 
of this field. More explicitly: at the start of the assumed MRI, 
the magnetic field has the value B,\ and then it exponentially 
grows to reach its maximum value Bm, followed by a rapid ex¬ 
ponential decrease reaching the final value Ba . The exponential 
growth and decay rates of the magnetic field are assumed to oc¬ 
cur with the same characteristic time. In a more general treat¬ 
ment, we can assume that the characteristic times of the field 
increasing is different from that of the field decaying. However, 
both times must be of the order of Alfven times of ms. A simi¬ 
lar field behaviour but for an electric field has been discussed in 
electromagnetism for the decay of the electric dipole moment 
( Schantz , 1995). Expectably, after this abrupt birth field de¬ 
cay, there will be a subsequent field decay caused by Ohmic 
diffusion and/or other resistive processes, which occur on time 
scales much larger than those of the birth MRI. 

An exponential growth occurs when the growth rate of the 
function is proportional to the current value of this function: 
/(f) oc /(f), where the overdot means time differentiation. Anal¬ 
ogously, an exponential decay occurs when /(f) oc -/(f). For 
the case of a magnetic moment fj(t) the positive and neg¬ 
ative growth rates can be described by the relation /i(f) oc 
-yu(f) tanh(f/r fl ), where the hyperbolic tangent function has 
been introduced to describe the positive and negative growth 
rates. The time t„ is the associated characteristic time. Accord¬ 


ingly, the behaviour of the magnetic moment of a newly-born 
neutron star during a birth MRI is assumed to be described by 
the equation fi(t) = -(l/r a )ju(f) tanh(f/r a ), whose solution reads 

sech (f/rj, (2) 

where //\i is the maximum value of the magnetic dipole mo¬ 
ment satisfying /j( 0) = fi\i- For a neutron star of radius R, the 
magnetic moment n is related to the magnetic field B by means 
of /r(f) = B(t)R } ( Jackson , 1998), which can be used together 
with Eq. (2) to obtain the expected magnetic field law 

B(t) = B m sech (f/T a ), (3) 

where Bm = B( 0) is the maximum value of the magnetic field. 
The behaviour of the magnetic field in the assumed birth MRI 
is qualitatively shown in Fig. 1. Equation (2) implies the non¬ 
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Figure 1: Qualitative behaviour of the magnetic field of a newly-bom neutron 
star during the assumed MRI. The times — ta/2 = -3 ms and t\/2 = 3 ms cor¬ 
respond to the beginning and end of this birth MRI. The value of the magnetic 
field at the beginning and end of the MRI is B\ = 10 13 G. The maximum value 
of the magnetic field =3 x 10 16 G occurs when B(0) = Bm • 


linear equation jj 2 - 4(ft 4 /iu 2 - fr /t 2 ) - n 2 /t 4 = 0, which com¬ 
bines with y u(t)-B M R 3 sech (f/r fl ) to yield 
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Using the Larmor formula P ra d(0 = 2/i(f) 2 /(3c 3 ) and Eq. (4), 
we obtain the instantaneous radiated power by a newly-born 
neutron star during the birth MRI, 


2 B 2 R 6 

p -=W s “ h 




2 -|2 

1 . 


(5) 


The change of kinetic energy of differential rotation can be writ¬ 
ten as d(adID. 2 /2)/dt = ( 4aan 2 MR 2 /5)d(l/P 2 )/dt. Using this 
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relation together with Eqs. (1) and (5) we obtain the instanta¬ 
neous energy conversion occurring in the birth MRI: 


values of the periods P a and Pj,. In fact, Eq. (10) with a d = 
0.1, M- 1.4M 0 and R- 10 km imply 
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Let us consider the time r A during which the MRI evolves. We 
will call r A the MRI time. We adopt the conditions B(±t a /2) = 
B A , where B A is the value of the magnetic field at the begin¬ 
ning and end of the birth MRI. These conditions and the law in 
Eq. (2) imply the result B a = Bm sech | t a /(2t„)] or equivalently 


T A =2T„asech(B A /B M ), (7) 


where asech denotes the inverse hyperbolic secant function. 

Integration of Eq. (6) over the time interval in which the as¬ 
sumed birth MRI evolves: -t a /2 < t < t a / 2, or equivalently, 
—T a asech (B A /B M ) < t < T a asech (B A /B M ), yields 
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where P a = P(—r A /2) and Pi, = P(+r A /2) are the values of 
the period at the beginning and the end of the birth MRI. The 
quantity v = v(t a /2) denotes the space velocity at the end of 
this instability, which is assumed to be the current value of the 
observed kick velocity. The initial condition v(-t a / 2) = 0 has 
been adopted, which is consistent with the idea that newly born 
stars acquire their space velocity during the birth MRI. 

If B M » B,\ then the radiation term in Eq. (8) reduces to 
AErad ~ (28/45)R 6 B^ i /(c 3 t 3 ) and thus Eq. (8) becomes 


P b= T 1 (ID 

AjPa 2 -(k\B 2 m +k 2 v 2 ) 

where k\ ~ .281 x ICE 30 cm/gr and k 2 ~ .633x 10~ 12 cm -2 . Since 
Pi, is real it follows that 

Pa< 1 = ■ (12) 

s jK l B 2 M +K 2 V 2 

This condition restricts the values of Bm and v. By assuming 
Bm = 3 x 10 16 G, we will apply Eqs. (11) and (12) to three dif¬ 
ferent velocities: (a) If v = 100 km/s then Eq. (12) implies 
P a < .0561 s. In particular, if P a = .015 s then Eq. (11) gives 
Pb ~ .0155 s; (b) If v = 500 km/s then P a < .0233 s. For exam¬ 
ple, if P a — .01 s then Pb ~ .0108 s; and (c) If v = 1000 km/s then 
P a < .0123 s. For instance, if P a = .005 s then Pb ~ .0055 s. If the 
initial period satisfies P a <- 015 s then the changes A P = Pb—P a 
(for the previous examples) are of the order of 10 4 s for the 
interval of space velocities 100 km/s < v < 1000 km/s. 

We can use Eq. (10) to obtain a formula for the kick velocity 
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where we have written v = v kick to emphasize that the kick ve¬ 
locity of neutron stars is originated by the birth MRI quantities: 
ad,Byi,Pa and /V This result explains the absence of the sus¬ 
pected correlation between the current magnetic field and space 
velocity, the so-called “v-B correlation” (Lorimer et ah, 1995, 
1997; Cordes & Chernoff , 1998). Let us rewrite Eq. (13) as 


Vkick yf^ikot} V jrad] ’ (^) 

where V| ( | mt | is the component of the kick velocity originated by 
the change of kinetic energy of differential rotation 
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V[drot] = -t 
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We expect that the time r A is of a few ms for magnetic fields of 
the order of 10 16 G. This expectation restricts the values of the 
time r fl and those of the relation B a /Bm- In particular, if r fl is 
taken to be r fl ~ 10 R/c then Eq. (7) implies a MRI time t a in 
the interval 3 ms < r A < 6 ms when Bm = 3 x 10 I6 G and B A 
lies in the interval 10 13 G<B A < 10 15 G. The choice r„ ~ 10 R/c 
is then consistent with the model proposed here. Insertion of 
T fl « 10 R/c into Eq. (9) yields 

7R 3 B^ M v 2 _4a d n 2 MR 2 j 1 1\ 
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and V[ la d] is the component of the kick velocity originated by 
the change of radiation energy 
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We can infer the kick force Fkick(0 producing v k i c k(0- Time- 
integration of Eq. (6) from -t a /2 to the time t leads to 


Vkick (0 = ^[droqM-V^M, (17) 

where 


Let us emphasize that Eq. (10) is valid when the magnetic field 
Bm is given by Eq. (3) and the conditions Bm » B A and T a ~ 
10 R/c are fulfilled. Equation (10) can be used to constrain the 
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in which k = (d>a d 7r 2 R 2 )/5, P a = P(-T\/2), and 
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If Eq. (18) is evaluated in the time t = t a /2 and we write 
P b = P( Ta /2) then we obtain Eq. (15). Analogously, if Eq. (19) 
with P la d given by Eq. (5) is evaluated in t — ta/2 and the con¬ 
ditions Bm » B,\ and r a ~ 10 R/c are fulfilled then we obtain 
Eq. (16). Using Eqs. (18) and (19) we can derive the corre¬ 
sponding accelerations. They can be expressed as 
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Time differentiation of Eq. (17) and the use of Eqs. (20) gives 
the kick (total) acceleration: a [kick] = v [d rot]a[drot]/v [k ick] - 
v [rad]a[rad]/ v [kickp which can alternatively be expressed as 
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We can evaluate a^y at t = 0 to get an idea of the order of its 
magnitude. According to Eq. (5), the power radiated at t = 0 is 
P r ad(0) « .1874 x 10 52 ergs for B u = 3 x 10 16 G, R = 10 km and 
T a = 3.33 x l(U 4 s. From Eq. (17) it follows that v [k i ck ](0) ~ 128 
km/s if we assume P a = .019 s, P( 0) = .01986 s and M— \AM 0 . 
We write .158 x 10 13 cm 2 (for a d = 0.1 and R- 10 km) and 
P( 0) «6.62 (calculated by P~2P/t a ). Using all of these values 
in Eq. (21) we get: a[ k ; ck ](0)«.53 x 10 8 g ( Heras , 2013). 

The forces associated with accelerations in Eqs. (20) read 
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By making use of F [kick] = (v [dl - ot] F [drot] -v [rad] F [rad] )/v [kick] and 
Eqs. (22) we can obtain the kick (total) force 
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It should be noted that Eq. (23) can naturally be interpreted 
in the context of the radiation reaction theory. If Ma = 
F|kick|> Text — kMP/(P V|kick|)' and 7v ad — Pi-ad/^[ k i ck ] then 
we obtain the well-known equation for the radiation reaction 
( Jackson , 1998): Ma = + ^rad- 

The velocity V[ ra dj in Eq. (16) can be related to the Alfven ve¬ 
locity V[Aif] = B/ A j4np. For a sphere of mass M and radius R the 
Alfven velocity reads V[ai f\=B/ -^2M/R 2 , which combines with 
Eq. (16) (making B = B M ) to yield V[ ra( i] = V2lv [A if]/75. This 
implies V[ A iq ~ 1 6.37v[ rad ]. If V[ ra d] ~ 200 km/s for B M = 3xl0 l6 G, 
M— 1.4M 0 and R- 10 km then V[ A if] ~3.27 x 10 8 cm/s, and this 
can be combined with \j A ir| =2R/T A (considering R = 10 km) to 
get the expected Alfven time T A ~6 ms, which emphasizes the 
consistency of the model proposed here. 

Let us apply Eqs. (14)-(16) to two representative neutron 
stars: The Crab pulsar B0531+21 and the magnetar J1809- 
1943. The canonical values M = 1.4M G and R- 10 km will be 


assumed. Unless otherwise specified, all data is taken from the 
ATNF Pulsar Catalogue (Manchester et ah, 2005). The space 
velocity v is obtained from the transverse velocity Vj_ by means 
of v=s V3/2 v k (Hobbs et ah, 2005; Lyne & Lorimer , 1994). 

Consider first the Crab pulsar with its current transverse ve¬ 
locity v ± = 141 km/s, which implies the space velocity v ~ 172 
km/s. It has been suggested that this pulsar was born with a 
period of 19 ms (Lyne et ah, 1993). If during a birth MRI in 
the Crab pulsar there was a very small increase in the period, 
for example, from P a = .019 s to P b = .02072 s and an abrupt 
change of magnetic field starting with the value 10 13 G, reach¬ 
ing the maximum value 3 x 10 16 G and ending with the value 
10 13 G, then Eq. (14) with a d = 0.1 yields Vkj ck 172 km/s 
and Eqs. (15) and (16) respectively give V[ drot ] ~ 264 km/s and 
v[r a d] ~200 km/s. Notice that the difference between the periods 
P a and P b is small: A/ J ~ .001 s. The current period of the Crab 
pulsar is P - .033 s. Crab’s period has increased about 13 ms 
since the occurrence of the assumed birth MRI. 

Consider now the magnetar J1809-1943, whose current pe¬ 
riod and magnetic field are P = 5.54 s and B - 2.1 x 10 14 G. 
This magnetar has v ± = 229 km/s, which implies v « 278 km/s. 
If this magnetar experienced a birth MRI in which there was a 
very small increase of its period, for example, from the value 
P a — .010 s to P b = .01039 s and an abrupt change of mag¬ 
netic field starting with the value 10 15 G, reaching the maximum 
value 3x 10 16 G and ending with the value 10 15 G, then Eq. (14) 
with a d = 0.1 yields V]d ck ~ 278 km/s and Eqs. (15) and (16) 
respectively give V[d ro t] ~ 342 km/s and V[ ra d] ~ 200 km/s. 

It is interesting to note that Eq. (11) allows the final period 
P b to have values considerably larger than those previously as¬ 
sumed. This possibility could explain why some neutron stars 
exhibit current periods in the scale of seconds (e.g., magnetars). 
According to Eq. (11) the period P b may be of the same order of 
the current period P. Consider, for example, the magnetar SGR 
1806-20 whose current period and magnetic field are P-1.5A s 
and B = 2x 10 15 G. This magnetar has v ± « 350 km/s (Olausen 
& Kaspi , 2014) which implies v ~ 428 km/s. If we assume 
the birth MRI values Z? A = 4 x 10 15 G and P a = .0265996 s then 
Eq. (11) implies P b = 6.01 s. This result suggests that shortly 
after the evolving of the assumed birth MRI, this magnetar ac¬ 
quired a period of the order of seconds. However, this conclu¬ 
sion faces the problem of how to physically explain the birth 
MRI change A P-P b - P„ of the order of seconds. 

In summary, any convincing explanation for the observed 
large space velocity of neutron stars should be traced to physical 
processes occurring during their birth. Accordingly, the space 
velocity should be related with birth values of physical prop¬ 
erties of these stars and not with current values of them. Here 
we have suggested that a MRI produced at the end of the birth 
process of neutron stars can be responsible for their large space 
velocities. We have shown how a rapid birth MRI conversion of 
kinetic energy of differential rotation into radiation energy and 
kinetic energy of translation occurring at the Alfven time of a 
few ms yields an observed interval for neutron star velocities, 
ranging from several hundreds to a few thousands km/s. 
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